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Abstract: Binuclear polypyridine ruthe- 
nium compounds have been shown to 
slowly intercalate into DNA, following 
a fast initial binding on the DNA sur- 
face. For these compounds, intercala- 
tion requires threading of a bulky sub- 
stituent, containing one Ru^\ through 
the DNA base-pair stack, and the ac- 
companying DNA duplex distortions 
are much more severe than with inter- 
calation of mononuclear compounds. 
Structural understanding of the process 
of intercalation may greatly gain from 
a characterisation of the initial interac- 
tions between binuclear Ru" com- 
pounds and DNA. We report a struc- 
tural NMR study on the binuclear Ru" 
intercalator A,A-B (A,A-[|Li-bidppz- 
(bipy)4Ru2]^+; bidppz = ll,ll'-bis(di- 
pyrido[3,2-a:2',3'-c]phenazinyl, bipy = 
2,2'-bipyridine) mixed with the palin- 
dromic DNA [d(CGCGAATTCGCG)]2. 



Threading of A,A-B depends on the 
presence and length of AT stretches in 
the DNA. Therefore, the latter was se- 
lected to promote initial binding, but 
due to the short stretch of AT base 
pairs, final intercalation is prevented. 
Structural calculations provide a model 
for the interaction: A,A-B is trapped in 
a well-defined surface-bound state con- 
sisting of an eccentric minor-groove 
binding. Most of the interaction enthal- 
py originates from electrostatic and 
van der Waals contacts, whereas inter- 
molecular hydrogen bonds may help to 
define a unique position of A,A-B. Mo- 
lecular dynamics simulations show that 
this minor-groove binding mode is 
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stable on a nanosecond scale. To the 
best of our knowledge, this is the first 
structural study by NMR spectroscopy 
on a binuclear Ru compound bound to 
DNA. In the calculated structure, one 
of the positively charged Ru^^ moieties 
is near the central AATT region ; this is 
favourable in view of potential interca- 
lation as observed by optical methods 
for DNA with longer AT stretches. Cir- 
cular dichroism (CD) spectroscopy sug- 
gests that a similar binding geometry is 
formed in mixtures of A,A-B with nat- 
ural calf thymus DNA. The present 
minor-groove binding mode is pro- 
posed to represent the initial surface 
interactions of binuclear Ru" com- 
pounds prior to intercalation into AT- 
rich DNA. 
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Introduction 

Structural investigations of interactions between DNA and 
small molecules are of paramount importance for develop- 
ing and improving chemotherapeutic agents that recognise a 
specific binding site.^^"^^ In this context, the unique proper- 
ties of octahedral-coordinated chiral (left-handed (A) or 
right-handed (A)) ruthenium(II) polypyridyl compounds 
have for many years captured great interest regarding DNA 
binding as optical probes, photo-reagents and inhibitors of 
DNA-related cellular processes. ^"^"^^ An example is the well- 
known monomeric prototype "light-switch" compound [Ru- 
(phen)2dppz]^^ (the monomer of compound 1 shown in 
Scheme 1; phen = l,10-phenanthroline; dppz = dipyrido[3,2- 
fl:2,3-c]phenazine), which intercalates into DNA by inserting 
its dppz moiety between the base pairs of DNA resulting in 
a thousand-fold enhancement of the luminescence bright- 
ness. Recent crystal and NMR structure studies of mono- 
meric Ru" compounds interacting with oligonucleotide du- 
plexes suggest that intercalation is sequence dependent, and 
that additional binding sites on the DNA may be exploit- 
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Scheme 1. Chemical structure of [|i-bidppz(phen)4Ru2]^^ (bidppz = 11,11 '- 
bis(dipyrido[3,2-a:2',3'-c]phenazinyl), a semi-rigid binuclear derivative of 
the "light switch " monomer [Ru(phen)2dppz]^''". 

Larger as well as more positively charged binuclear com- 
pounds have been observed to be more efficient in sensing 
larger scale sequence-specific textural properties of the nu- 
cleic acid structure The binuclear Ru" compound 
bidppz(phen)4Ru2]'^^ (1) (Scheme 1) and its derivatives have 
been shown to intercalate by inserting the bridging bidppz 
ligand between base pairs, leaving one Ru(phen)2 moiety in 
each groove. However, the dumbbell shape of these dimers 
(Figure 1) requires a very different mode of intercalation 




Figure 1. Structure of A,A-[|i,-bidppz(bipy)4Ru2]'^'^ (A,A-B). Colour 
coding: carbon atoms black, nitrogen atoms grey and hydrogen atoms 
white ; Ru^^ ions are shown as spheres. 



compared to the monomer: instead of shding a planar struc- 
ture, such as dppz of a monomeric Ru" compound, between 
two base pairs, the bulky substituent surrounding a Ru" ion 
has to be threaded through the base-pair stack; this requires 
severe DNA distortions and probably also transient opening 
of one base pair.^^^^ Just like the monomer, intercalation of 
the bidppz ring system of binuclear Ru" compounds is ac- 
companied by a very large luminescence increase. ^^^^ Upon 
interaction with calf thymus DNA (ct-DNA) the A,A-enan- 
tiomer of 1 has been shown by linear dichroism (LD) to re- 
organise from an initial surface-bound geometry into an in- 
tercalative binding mode in a very slow process (two weeks 
at room temper ature).^^^^ 

Ru" compounds such as 1 have been shown to selectively 
thread-intercalate into long AT-rich regions, ^^^^ the 2,2'-bi- 
pyridine analogue of compound 1, A,A-[|Li-bidppz- 
(bipy)4Ru2]'^^ (A,A-B; Figure 1) being the most discrimina- 
tive. ^^^^ Although threading into [poly(dAdT)]2 occurs within 
minutes, this interaction rate is about 2500-fold slower with 
natural DNA such as calf thymus DNA. A certain length of 
the stretch of AT base pairs is required for fast intercala- 
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tion: a dramatic drop in the threading rate was observed 
when the number of consecutive base pairs (bp) with alter- 
nating AT sequence was decreased from 14 to 10 bp, and no 
threading seems to occur for a stretch of six AT base 
pairs, f^^^ Another significant observation concerns the inter- 
calation rate. This is also in the best cases slow, varying at 
room temperature from a ty2 of 3.4 min for [poly(dAdT)]2 
and weeks for calf thymus DNA.^^^^ Thus, the intercalation 
process must occur in two steps. Diffusion guided by the 
electrostatic attraction by the DNA phosphate groups re- 
sults in the rapid formation of a complex between A,A-B 
and the DNA, where A,A-B is bound to the surface of the 
DNA, perhaps near an AT-rich region due to the local in- 
crease of the electronegative potential surrounding these re- 
gions. ^^"^^ Provided a suitable thermodynamic interaction 
force is present, this is then followed by a slow re-orienta- 
tion eventually yielding an intercalating state. Although this 
latter state has been detected by optical methods, ^^^^ little is 
known about the process of intercalation and in particular 
about the intermediate surface bound state. LD studies esti- 
mate that the angle between the long axis of A,A-B and the 
DNA helix axis is about 65°.^^^^ However, no 3D structure of 
a binuclear Ru" compound bound to DNA has been report- 
ed so far. 

By characterising the intermediate, surface-bound geome- 
try, the following questions can be addressed: in which 
groove does A,A-B initially bind and therefore, from which 
groove does it thread-intercalate into DNA; how do the 
bulky "propellers" of the rather rigid A,A-B (Figure 1) 
orient themselves with respect to the DNA, and what inter- 
actions stabihse this initial complex. Here, an NMR structur- 
al study based on NOE observations is combined with MD 
(molecular dynamics) simulations to test the (short-term) 
stability of the resulting structures. The selection of a DNA 
duplex was based on the following requirements: it should 
contain an AT-rich region that attracts A,A-B, but the se- 
quence should not allow any intercalation, but rather trap a 
surface-bound intermediate state. The well-studied palin- 
dromic 12-mer [d(CGCGAATTCGCG)]2, originally con- 
structed by Dickerson et al. fulfils these demands, and 
offers several CG base pairs for stabilisation of the ends. 
Note however, that the short AT region does not allow both 
Ru^"^ centres of A,A-B to simultaneously approach this cen- 
tral stretch. Thus, the additional question arises as to how 
the two Ru^"^ centres adjust to the presence of only four A- 
T base pairs. 

Results 

Optical spectroscopy: Circular dichroism (CD) and lumines- 
cence spectra were compared for various DNA sequences 
mixed with A,A-B (Figure 2). At room temperature, A,A-B 
thread-intercalates into [poly(dAdT)]2 within less than 
5 min, but only very slowly into calf thymus DNA (ct- 
DNA). ^^^'^^'^^^ The CD spectra of both 
[d(CGCGAATTCGCG)]2 (Figure 2, green curve) and ct- 
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Figure 2. Circular dichroism spectra of mixtures with 50 |im A,A-B in 
20 mM phosphate buffer (pH 6.5), 100 mM NaCl water solution: 
[d(CGCGAATTCGCG)]2 (green curve), ct-DNA (red curve, mostly cov- 
ered by the green curve) and [poly(dAdT)]2 (blue curve). The [nucleoti- 
des]/[A,A-B] ratio is 24:1 for all mixtures. Inset: Luminescence spectra of 
the corresponding mixtures with the same colour coding, excitation wave- 
length: /l = 410nm. All spectra were recorded about half an hour after 
mixing at 25°. 



DNA immediately after mixing (Figure 2, red curve) with 
A,A-B are strikingly similar, but differ strongly from the CD 
spectrum of the [poly(dAdT)]2/A,A-B mixture (Figure 2, 
blue curve) J^^^ These pronounced spectral differences sug- 
gest that AA-B binds to [d(CGCGAATTCGCG)]2 in a sim- 
ilar way as the initial binding to ct-DNA J^^^ Examination of 
the same samples with fluorescence spectroscopy shows a 
six times lower magnitude of luminescence intensity for 
both ct-DNA and 

[d(CGCGAATTCGCG)]2 when 
compared to [poly(dAdT)]2, the 
reduced luminescence reflecting 
access to solvent water for ex- 
ternally bound complex,^^^^ thus 
supporting the conclusion from 
the CD data. 



with A,A-B binding to DNA. With further titration, the 
excess of A,A-B yielded strong line broadening, preventing 
observations of any other changes in the spectrum (Fig- 
ure S2 in the Supporting Information). Therefore, a 1:1 ratio 
of A,A-B and [d(CGCGAATTCGCG)]2 was selected for all 
2D NMR experiments. It proved necessary to add A,A-B to 
the DNA solution; the reverse order resulted already at 
very low concentrations of DNA in precipitation and disap- 
pearance of signals. 

Assignment of free DNA and free A,A-B : NMR reso- 
nances of free [d(CGCGAATTCGCG)]2 were assigned by 
using standard techniques.^^^'^^^ Five imino proton resonances 
were observed, indicating that the oligonucleotide forms a 
stable double helix at the given conditions, with fraying indi- 
cations only for the first and last base pairs, observations 
that are in agreement with previous studies. ^^^'^^^ A B-type 
conformation is confirmed by the position of the two most 
downfield-shifted resonances, the imino protons for the thy- 
mine bases ((5 = 13.75 and 13.61 ppm), followed by those of 
the guanine bases ((5 = 13.05, 12.88 and 12.68 ppm). Com- 
parison of the chemical shifts with literature values^^^^ re- 
veals only small deviations (Table SI in the Supporting In- 
formation). 

Free A,A-B in 20 mM phosphate buffer solution at room 
temperature was assigned by analysing TOCSY and 
NOESY spectra. One three-spin system and four four-spin 
systems were detected. Figure 3 (lower panel) shows the res- 
onance assignment in a ID spectrum (chemical shifts are 



Conditions for NMR measure- 
ments: Comparison of NOESY 
spectra for samples with a 1:1 
mixture of 
[d(CGCGAATTCGCG)]2 and 
A,A-B, in 100 mM NaCl solution 
or in 100 mM NaCl solution 
with 20 mM sodium phosphate 
buffer (pH 6.5), showed that 
without buffer several imino 
proton peaks for the GC base 
pairs were missing and several 
cross peaks became weaker 
than the corresponding ones in 

buffered solution (Figure SI in the Supporting Information). 
Therefore, 20 mM sodium phosphate buffer with 100 mM 
NaCl, pH 6.5, was utilised in all samples. When titrating 
A,A-B into a DNA solution starting from a ratio of 0:1 to a 
final ratio of 1:1, new resonances were appearing consistent 




\1 8.0 7.8 
Chemical shift (ppm) 

Figure 3. NMR assignments of 1 mM A,A-B in 20 mM phosphate, pH 6.5, 100 mM NaCl water solution. Res- 
onance assignments refer to positional symbols defined in the chemical structure of A,A-B with positional 
symbols at the carbon atoms at the top. Initial numbers (with or without primes) indicate the position within 
the rings, which are identified by capital letters. "A" and "B" are the axial and equatorial (with respect to the 
dppz plane) pyridine rings of the bipyridine ligands, respectively; "C" denotes the pyridine and "D" denotes 
the benzene rings of the dppz hgand. Primed and non-primed numbers discriminate the two bipyridine systems 
attached to the ruthenium atoms. Lower case letters define the two equivalent monomer units. 



listed in Table S2 in the Supporting Information). Due to 
the overall two-fold symmetry, both monomer units, labelled 
"a" and "b" in Figure 3, are identical. The two axial (A, A') 
and the two equatorial (B, B') pyridine rings show slightly 
different chemical shifts due to the breaking of the mono- 
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mer two-fold symmetry by the pivot bond, however, the two 
pyridine rings (C and C) of the bidppz hgand appear identi- 
cal. The extraordinary downfield-shifted peak (6 = 
9.86 ppm) was assigned to H4C/H4'C because of their close 
distance to the electron-withdrawing nitrogen atoms of the 
phenazine part of the bidppz ligand. 

Assignments of bound DNA and bound A,A-B : Upon inter- 
action between A,A-B and [d(CGCGAATTCGCG)]2, the 
DNA signals are split into three sets, each of which can be 
assigned to an individual, complete or nearly complete 
DNA strand. The strongest signals can be assigned to a sym- 
metric duplex, very similar to the one for free DNA, and a 
basically complete assignment was achieved. Besides this 
symmetric DNA duplex, somewhat weaker signals charac- 
terise two additional strands with almost complete assign- 
ment (in the following referred to as a strand and (3 strand), 
which together form an asymmetric DNA duplex based on 
interstrand NOEs. Sequential walks for these two strands 
are shown in Figure 4: green for the a strand and orange for 
the P strand (sequential and interstrand NOEs are identified 
in the structure of Figure S3 in the Supporting Information). 
On average, the (3 strand shows stronger chemical shift dif- 
ferences from free DNA than the a strand (Figure S4 in the 
Supporting Information). In Figure 5, the region with the 
imino proton signals is shown for all three DNA strands, by 
using the same colour coding as in Figure 4 for the asym- 
metric strands. Although these NOE networks are incom- 



3-8T3-7T 1-lOG 
3-8T 3-7T 3-7T 3-8T 1-2G 1-lOG 



1-4G 
1-4G 




41-3C 



13.6 13.4 13.0 12.8 12.6 
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Figure 5. Imino region of the NOESY spectrum shown in Figure 4. Cross- 
peak assignments for all three observed DNA strands are given with 
numbers indicating the position in the base rings, followed by the nucleo- 
tide number and type: symmetric DNA (blue), asymmetric DNA with 
green (a strand) and orange ((3 strand). Vertical and horizontal lines con- 
nect peaks to the assignment texts. 



8-5A 8-6A 8-6A 8-lOG 8^ 8-4G 



6-9C 6-8T 
6-9C 6-8T 6-3C 6-3C 6-7T 



r-5A r-5A r-7T l'-9C l'-4G 

r-8T r-8T r-6A v-n r-3C r-4G 




73 TA 6.9 63 

Chemical shift (ppm) 

Figure 4. Selected regions of the NOESY spectrum of a 1:1 mixture of 1 mM A,A-B with 
[d(CGCGAATTCGCG)]2, showing connections between base and sugar protons. Sequential walks are pre- 
sented for the a strand (green) and the (3 strand (orange). Individual cross peaks are identified on the left and 
top borders with numbers indicating the positions in the sugar or base rings, followed by the nucleotide 
number and type. 



plete, lacking connections to 
the terminal nucleotides, their 
observation clearly excludes the 
possibility of intercalation any- 
where between base pairs 2 and 
9 (a complete hst of the chemi- 
cal shifts for all DNA strands as 
well as further assignment in- 
formation is given in Table SI 
in the Supporting Information, 
and chemical shift differences 
of the a strand and the (3 strand 
are shown in Figure S4 in the 
Supporting Information). 

In the NMR spectra of the 
complex between A,A-B and 
[d(CGCGAATTCGCG)]2, four 
resonances of A,A-B are found 
with chemical shifts d>9 ppm, 
that is, outside the spectral re- 
gions occupied by DNA reso- 
nances (Figure S7 in the Sup- 
porting Information). Two reso- 
nances exhibit intermolecular 
NOEs (see below) to the DNA: 
one appears at 6 = 935 and the 
other at 9.14 ppm. Based on the 
patterns of intramolecular 
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NOESY and TOCSY peaks and chemical shift differences 
to the resonances of free A,A-B (for details see Section 3 in 
the Supporting Information), these four resonances can only 
be assigned to the H4Ca, H4'Ca, H4Cb and H4'Cb protons; 
these are expected to have individual resonances after sym- 
metry breaking. Unambiguous individual assignments of the 
four atoms to these resonances are not possible at this stage, 
and have to await structural calculations. 

Intermolecular NOESY cross peaks and CYANA calcula- 
tions: The assignment of the DNA resonances for intermo- 
lecular NOEs between A,A-B and the asymmetric DNA 
duplex is illustrated in Figure 6 : Each NOE is connected to 
the network of intramolecular DNA NOEs used for the 
DNA assignments (Figure 4). As discussed above, only the 
H4Ca/b and H4'Ca/b protons remain as candidates for the 
assignment of the two resonances at 6 = 935 and 9.14 ppm 
in the spectra of the mixture between A,A-B and 
[d(CGCGAATTCGCG)]2. From all possible assignment 
combinations (four available protons for each resonance 
yield in total sixteen possibilities), the four with the same 
proton assigned to both resonances can immediately be ex- 
cluded. The remaining twelve combinations consist of identi- 
cal pairs due to the symmetry of A,A-B. For the remaining 
six combinations, CYANA calculations were performed. For 
one assignment combination, the CYANA calculations show 
no violations of the distance restraints or of any van der 
Waals limits exceeding 0.1 A combined with CYANA target 
functions as low as 0.13. The other five assignment combina- 
tions resulted in violations of both distance and van der 
Waals restraints of at least 0.35 A, and target functions ex- 
ceeding 0.43 (Table S3 in the Supporting Information). Thus, 
the assignment with smallest violations and lowest target 
function was accepted, with the following three NOEs: 
proton 4Ca of A,A-B to the HI' protons of both G4 and A5 
(|3 strand), and proton 4'Cb of A,A-B to the HI' proton of 
both G2 (P strand). The resulting five best structures calcu- 
lated with CYANA (see the Experimental Section) by using 
these three NOEs show that A,A-B is aligned with the 
minor groove with the two monomers in an anti conforma- 
tion with a torsion angle of 33°. 

Model structure refinement with molecular dynamics: The 

resulting CYANA model was subjected to exhaustive MD 
simulations with the aim of further refining the structure, to 
evaluate its (short-term) stability and to confirm the pres- 
ence of a stable minor-groove binding mode of A,A-B that 
satisfies, on a time average, the intermolecular NOEs. The 
best of the five resulting CYANA structures with the accept- 
ed NOE assignment served as starting structure for ten fully 
unrestrained MD runs, each 5 ns long, after equilibration 
and REMD (see the Experimental Section). Consistency 
with the NOE restraints was tested by calculating for each 
observed intermolecular NOE an average distance (d^oE) 
between the two protons according to Equation (1) (see the 
Experimental Section). Table 1 provides comparisons be- 
tween the NOE distance limits and the corresponding aver- 
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4Ca 8-5A l'-5A l'-4G 




(b) 

Chemical shift (ppm) 

Figure 6. Spectral regions from the NOESY spectra shown in Figures 4 
and 5. a) NOEs between the A,A-B proton 4Ca (see Figure 3) at d = 
9.14 ppm to HI' of 4G and HI' of 5A. b) NOE between the A,A-B 
proton 4'Cb at (5 = 9.35 ppm to HI' of 2G. All NOEs are to the (3 strand. 
To confirm the DNA assignments, these intermolecular NOEs are con- 
nected to the intramolecular network of the NOEs of the DNA. The 
cross peaks are identified on the left and top borders. 

age distances d^o^ for each MD run. Complete consistency 
is confirmed for five runs (4-7 and 10). For run 2, the maxi- 
mal violation does not exceed 0.1 A. Remarkably, for all ten 
runs the two shorter restraints with limits of 3.5 A are fulfil- 
led on a time average (with the above small exception), and 
only the longer restraint is sometimes violated. This latter 
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Table 1. Assigned NOEs used for the structural calculations. DNA signals are all from 
the [3 strand. Distances <iNOE according to Equation (1) are given for the ten unre- 
strained MD runs of 5 ns each. The five runs that are strictly consistent with the 
NOEs are at positions 4-7 and 10. 



NOE 




A,A-B 


Limit^^] 






d^oE in 10 MD runs [A] 






1 


G4 


4CA 


5.0 


5.4 


4.2 


5.7 4.9 5.0 4.6 4.6 5.5 


5.7 


4.6 


2 


A5 


4CA 


3.5 


2.7 


2.8 


2.9 2.7 2.8 2.7 2.8 3.1 


2.8 


2.7 


3 


G2 


4'CB 


3.5 


3.2 


3.6 


3.4 3.0 3.1 3.0 2.9 3.5 


2.4 


2.9 



[a] HI' protons from the (3 strand of the DNA. [b] Distance restraints from NOE ob- 
servations in [A] . 



tance of 4.4 A from the centre of the A,A-B ring 
with atoms *'Cb (Figure 3), and the angle between 
the connection of this proton to the ring centre and 
the ring plane is 7°. This position near the ring 
plane and close to the ring explains the shift change 
as a ring current effect. The H8 protons of A5 and 
A6 as well as the H6 proton of T7 on the other 
hand are in the major groove and at distances of 
6.7, 7.2 and 9.0 A of the closest A,A-B ring; the cor- 



violations result from either a small displacement of the 
A,A-B along the minor groove or small changes in the DNA 
structure. In one simulation (run 3), structural divergence in- 
dications become apparent towards the end of the run. 

3D structure of the DNA-A,A-B complex: The fact that at 
least five runs are stable over 5 ns of unrestrained MD 
strongly supports the existence of a stable complex with 
A,A-B bound to the minor groove of the DNA in a manner 
consistent with all NOE restraints. For structural compari- 
sons, averaged structures were determined from the snap- 
shots of each MD. The relative positions of A,A-B in the dif- 
ferent MD trajectories were then compared as follows: the 
average structures of each MD run were superimposed for 
minimal root-mean-square deviation (RMSD) of the most 
stable parts of the binding half of the DNA, that is, base 
pairs 3-6 (all heavy atoms). Following this superposition the 
average displacement of A,A-B heavy atoms was (0.6 ± 
0.3) A. Visual inspection shows that most variations affect 
the outer bipyridine rings. This well-defined position of 
A,A-B in the minor groove is supported by the identification 
of two hydrogen bonds (length < 2.6 A, angle < 25°) in all 
five structures connecting Na and N'b, of A,A-B (see 
Figure 1) with the amino groups of the guanine bases 4 and 
2, respectively. 

The average structure of run 7 yielded the smallest A,A-B 
displacements to the other average structures, and is there- 
fore a good representative of all MD runs. Its A,A-B heavy 
atom displacements and the variations of the distances cor- 
responding to the NOEs during the unrestrained 5 ns MD 
run is presented in Figure 7, demonstrating again the stabili- 
ty of the complex. Figure 8 shows the location of A,A-B in 
the minor groove of the DNA of this average structure. The 
angle between the main axis of the dimer and the DNA 
helix axis is consistent with the LD value of about 65°.^^^^ 
One Ru^^ moiety is near the central A ATT region. This 
binding geometry is secured by the hydrophobic alignment 
of the bridging bidppz ligands of the A,A-B with the fura- 
nose rings of the DNA backbone, as well as by electrostatic 
interactions of the positively charged peripheral hydrogen 
atom of the ligand with the negatively charged phosphate 
groups of the DNA backbone. 

The binding of A,A-B to the DNA also results in chemical 
shift changes of the protons of the DNA, the largest four ex- 
ceeding 0.35 ppm (Figure S4 in the Supporting Information). 
The HI' of G2 is in the structure of Figure 8 located at a dis- 
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Figure 7. Averaged structural variations versus simulation time for MD 
run 7. Only the last unrestrained 5 ns are shown. The bottom three 
panels show the time evolution for proton-proton distances correspond- 
ing to the measured NOEs: a), b) and c) correspond to NOE 1, 2 and 3 
in Table 1, respectively, d) Average displacement of heavy atoms of A,A- 
B between the snapshot structures and the mean structure of the corre- 
sponding trajectory interval. For these displacements, the DNA heavy 
atoms of base pairs 3-6 have been superimposed, followed by calculation 
of position differences for only the heavy atoms of A,A-B. 



responding angles are between 30 and 60°. This excludes a 
ring current effect. However, changes of approximately 
0.5 ppm for base protons in the major groove upon binding 
of a minor-groove ligand have been observed earlier and are 
probably due to conformational changes of the DNA.^^"^"^^^ 

The interaction enthalpy was estimated by summing inter- 
molecular van der Waals and Coulomb interactions, and by 
adding contributions from the two hydrogen bonds. With 
calculations as described in the Experimental Section, one 
obtains estimates of -(64.5±6.3) and -(83.2±0.33) kcal 
mol"^ for van der Waals and Coulomb interactions, respec- 
tively (averaged over the average structures of the five 
stable and NOE-consistent MD runs); the two hydrogen 
bonds add only a few kcalmol"^ to this interaction enthalpy 
(note that this enthalpy estimate does not include any free 
energy contribution caused by hydrophobic effects, which 
for this highly hydrophobic compound A,A-B is expected to 
further favour complex formation). Due to the asymmetric 
position of A,A-B with respect to the DNA centre, one of 
the monomers interacts with one end of the DNA. Howev- 
er, only (1.6 ±0.6)% of the total van der Waals and Cou- 
lomb interaction enthalpy — (147.7±6.3) kcalmol"^ involves 
the DNA end surface (i.e., the surface created by cutting a 
12-mer fragment from a long DNA). 
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Figure 8. Average DNA-A,A-B structure of the 5 ns productive, unrestrained trajectory from the same MD simulation as shown in Figure 7. A) Overall 
side view. B,C) Details of the short distances of three intermolecular NOEs shown in Table 1. Colour coding: green = strand a, orange = strand (3. Dotted 
lines indicate NOEs (black, numbered 1-3) and potential hydrogen bonds (orange, labelled "a" and "b": Na and N'b of A,A-B with the amino groups of 
the guanines bases 4 and 2, respectively. 



Discussion 

In the NMR spectra, most of the A,A-B resonances overlap 
with the protons of the DNA bases at 6 = 7-8 ppm. This re- 
sults in a very limited number of assignable intermolecular 
NOEs. However, because all observed intermolecular NOEs 
are from the bidppz ligand to the sugar protons of the 
DNA, the position of the bidppz moiety in the minor 
groove can be well determined. Because the flexibility of 
A,A-B is reduced to the single bond connecting the two 
monomers, the position of the bipyridine rings can be also 
well defined. 

The present structure determination by NMR spectrosco- 
py describes a 1:1 complex where A,A-B binds to the minor 
groove near one end of the asymmetric DNA duplex. The 
NMR data show in addition the presence of a symmetric 
DNA duplex (free DNA) in the NMR sample, but no addi- 
tional form of A,A-B could be detected (see also the Sup- 
porting Information), with the exception of a fast exchange 
process located at the terminal DNA base pairs of the free 
DNA (Section 5 in the Supporting Information). In particu- 
lar, intercalation can be excluded based on the near-com- 
plete sequential NOE connections for the asymmetric DNA 
duplex (Figure S3 in the Supporting Information), in agree- 
ment with the non-intercalative binding mode suggested by 



the CD and luminescence results. Interestingly, the NMR 
data exclude the presence of a significant concentration of a 
symmetrically bound A,A-B at the other end of the asym- 
metric DNA duplex (e.g., near the bottom in Figure 8 A), al- 
though no steric clash would hinder such a 2:1 complex. The 
complexity in our NMR spectra, with signals from symmet- 
ric and asymmetric DNA forms, precludes a precision of the 
resulting DNA structure that would for example allow de- 
termining conformational changes preventing binding at the 
other end. To our knowledge, this is the first atomic resolu- 
tion structure of a complex between a binuclear Ru" com- 
pound and a DNA fragment. As discussed below, we pro- 
pose that this (rapidly attained) surface-binding mode ap- 
proximates the starting point for the (slow) intercalation 
process observed with AT-rich DNA through optical meas- 
urements. ^^^^ 

Besides the obvious electrostatic attraction between A,A- 
B with four positive charges and the strongly negatively 
charged DNA, the observed minor-groove structure is also 
stabilised by burying the highly hydrophobic surface of A,A- 
B on the DNA surface, allowing two intermolecular hydro- 
gen bonds to be formed to the phenazine nitrogen atoms. 
However, although these may help to define a unique posi- 
tion of A,A-B in the minor groove, enthalpy calculations 
(estimated as the sum of van der Waals and electrostatic in- 



Chem. Eur. J. 2013, 19, 5401-5410 



© 2013 Wiley- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.chemeurj .org 



5407 



CHEMISTRY= 



A EUROPEAN JOURNAL 



B. Norden et al. 



teractions and hydrogen bonds as explained in the Experi- 
mental Section) indicate the overwhelming intermolecular 
attraction comes from the sequence-independent van der 
Waals, electrostatic and hydrophobic contributions. A simi- 
lar argument shows that in spite of the location of A,A-B 
near one end of the DNA, the surface created by cutting the 
12-mer DNA contributes negligibly (about 1.6%) to the 
total interaction enthalpy. The observed preference of Ru" 
compounds for AT-rich DNA region^^^^ may explain the po- 
sition of one A,A-B monomer with a central Ru ion near 
the AT base pairs, for which an increase of the electronega- 
tive potential was described. ^^"^^ This results in the position of 
the second monomer at the end of the DNA duplex. When 
the upper end of the B-DNA model shown in Figure 8 is ex- 
tended by four base pairs, a steric clash will result between 
the left wall of the minor groove and the upper left bipyri- 
dine ligand. A subsequent refinement suggests that this 
clash can be resolved by a slight increase of the width of the 
minor groove and a small tilt to the right wall of the minor 
groove of the upper part of the complex while still fulfilling 
the NOE constraints. Superposition of the upper six base 
pairs of the model in Figure 8 with the corresponding base 
pairs of the longer DNA results in a RMSD of 0.8 A and an 
average displacement of the Ru" compound of 1.7 A, leav- 
ing the position and orientation of the lower Ru" monomer 
nearly unperturbed near the AATT region. 

Previous linear dichroism and luminescence studies re- 
ported that binding to the DNA surface is the first step of 
the slow thread-intercalation of binuclear Ru" compounds. 
Threading intercalation is observed for various binuclear 
Ru" compounds including A,A-B, provided a sufficiently 
long stretch of AT base pairs is present. The stable minor- 
groove binding of A,A-B to [d(CGCGAATTCGCG)]2, 
which contains a central sequence of AT base pairs that 
does not allow intercalation to occur, is therefore likely to 
correspond to a trapped intermediate prior to intercalation. 

Binuclear Ru" compounds differ from the better studied 
monomeric compounds, such as A-[Ru(phen)2dppz]^+ and 
A-[Ru(bpy)2dppz]^+, mostly by the need to thread the bulky 
substituents surrounding one Ru" ion (Figure 1) through the 
base-pair stack. In contrast to simply sliding a planar ring, 
for example, dppz, between two base pair as needed in the 
monomeric case,^^'^^^ a major DNA distortion, likely com- 
bined with opening of a base pair, is required. Nonetheless, 
it is interesting to notice that also the monomeric com- 
pounds prefer to intercalate from the minor groove as illus- 
trated by two recent crystallographic studies. ^^'^^^ Besides dif- 
ferent intercalation modes (threading between two base 
pairs, filling the hole created by a mismatch, entering the 
DNA at the fraying ends), these crystal structures also 
reveal subtle dependencies on the DNA sequences, for ex- 
ample, favouring AT over TA. 

Binuclear Ru" compounds such as A,A-B differ from 
both small molecules such as netropsin and many others, ^^^^ 
as well as proteins^^^^ binding in the minor groove. Practical- 
ly all DNA binding molecules utilise positive charges to 
offset the negative charges of the DNA backbone and the 



negative electrostatic potential in the grooves. On the other 
hand, compounds such as A,A-B are highly rigid, whereas 
small molecules can adapt their conformation to the curva- 
ture and shape of the bottom of the minor groove. Similarly, 
the flexibility of the side chains of proteins allows them to 
contact specific groups of the DNA bases. The ability of 
both small molecules and proteins to form hydrogen bonds 
to several hydrophihc groups of the DNA bases provides a 
higher potential for sequence specific DNA binding. Due to 
the rigidity of A,A-B, the adaptability to form hydrogen 
bonds is strongly reduced compared to small molecules or 
proteins. A,A-B and proteins share the size characteristics 
that prevent full integration in the minor groove. Together 
with the predominantly hydrophobic character and the posi- 
tive charges of A,A-B, this explains the unique binding 
mode of A,A-B with partial-groove penetration of only the 
central bidppz moiety, as well as the attraction to AT-rich 
regions, and possibly the ability to thread-intercalate into 
suitable DNA. 



Conclusion 

In summary, in the 1:1 DNA/A,A-B complex we have shown 
that A,A-B binds to the minor groove of 
[d(CGCGAATTCGCG)]2 in an asymmetric manner. Gener- 
al long-range electrostatic interactions are the driving force 
for the initial attraction between the two molecules; local 
variations of the electrostatic potential may then direct the 
Ru" compound near the AT base pairs, defining the asym- 
metric location on the DNA. Based on estimation of binding 
enthalpies, a similar minor-groove binding can also be ex- 
pected for other DNA sequences, including the ones with 
sufficiently long AT stretches that allow subsequent thread- 
ing intercalation. The strikingly similar CD spectra obtained 
for the complex of A,A-B with [d(CGCGAATTCGCG)]2 
and for freshly prepared mixtures of A,A-B with ct-DNA 
(before the slow intercalation process has taken place) indi- 
cate that a similar binding mode is present with the natural 
calf thymus DNA/A,A-B complex. To the best of our knowl- 
edge, this is the first NMR structure on binuclear ruthenium 
A,A-B interaction with DNA in a surface-bound manner. 
Future work will focus on a DNA sequence that allows de- 
termining the structure of a surface-bound mode before ini- 
tiating threading and structural analysis of the intercalated 
mode. 



Experimental Section 

Materials: The oligomer [d(CGCGAATTCGCG)]2 was synthesised with 
standard phosphoramidite chemistry, purified by HPLC, lyophihsed and 
subsequently re-dissolved in 100 mM NaCl solution, in the presence of 
20 mM sodium phosphate (pH 6.5, in the following referred to as phos- 
phate buffer), providing stock solutions for NMR and optical spectrosco- 
py. Calf thymus DNA (ct-DNA) and [poly(dAdT)]2 were obtained from 
Sigma-Aldrich and used as received. DNA concentrations were deter- 
mined by absorbance measurements on a Gary 4B spectrophotometer, by 
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using molar extinctions £260 = 193070 cm^M^ (calculated by using Oli- 
go Analyzer software, available online at http://www.idtdna.com) for the 
oligomer duplex, £260 = 6600 cm^ for ct-DNA and £262 = 6600 cm^ 
for [poly(dAdT)]2. Compound A,A-B, as tetrachloride salt, was synthes- 
ised as described elsewhere. ^^^'^^^ Its concentration was determined by 
using £410 = 65000 cm^M ^ 

Spectroscopy: CD spectra were acquired on a Jasco J-810 spectropo- 
larimeter by using a 1 cm quartz cell. CD samples were prepared 50 |im 

A, A-B in 20 mM phosphate buffer (pH 6.5), 100 mM NaCl solution. The 
[nucleotides]/[A,A-B] ratios are 24:1 for all mixtures, and CD spectra 
were recorded immediately after mixing. Emission spectra of the same 
samples were measured shortly after the CD scans on a Varian Eclipse 
spectrofluorimeter with excitation wavelength /l = 410nm. For all NMR 
measurements DNA was annealed by heating to 50 °C, where the imino 
signals disappeared, followed by slow cooling to 25 °C (1 °Cmin^), where 
the imino signals reappeared. All samples were buffered with 20 mM 
sodium phosphate and the pH value was adjusted to 6.5. Spectra were re- 
corded in Shigemi tubes with 1 mM DNA concentration in a 9:1 mixture 
of H2O/D2O. For the samples containing mixtures of DNA and A,A-B, 
the latter was added after the annealing. All NMR spectra were recorded 
at 25 °C on a Varian Inova 800 MHz spectrometer by using Watergate 
solvent suppression. One-dimensional spectra were recorded with 16 K 
data points; two-dimensional NOESY and TOCSY spectra were acquired 
by using 512 data points in ti and 2048 data points in ?2, with a pulse se- 
quence repetition delay of 1.5 s. Both types of 2D spectra were recorded 
on pure DNA and on 1:1 mixture of the DNA with A,A-B, by using 
mixing times of 100 ms for NOESY and 80 ms for TOCSY experiments. 
All data were processed with NMRPipe^^^^ by zero-filling to 1024 points 
in and apodising with either a Gaussian or a shifted sine-bell function. 
XEASY^^^^ and CCPNmr^^^^ were used for analysis, assignment and peak 
integration of 2D spectra. ID spectra were analysed with Mnova.^^^^ 
Structure calculations: Different NOE inputs for the complex with A,A- 
B and [d(CGCGAATTCGCG)]2 were tested with CYANA^^^] calcula- 
tions by using rigid B-DNA strands, whereas for the A,A-B dimer a 
single torsion angle was defined for the bond connecting the two mono- 
mer units (Figure 2). The two DNA strands were connected by a flexible 
linker comprising 14 standard CYANA hnker units (these hnker units 
consist of "pseudo-atoms" that have zero van der Waals radius); similar- 
ly, the second DNA strand was connected to the A,A-B with a 19 units 
long linker. A standard B-DNA double helix structure was enforced by 
hydrogen-bond restraints for all twelve base pairs. In addition, observed 
NOEs were used to define upper distance limits between A,A-B and 
DNA, by using upper limits of 5 A for weak and 3.5 A for strong NOEs. 
The standard CYANA protocol^^"*^ was applied to calculate 50 structures, 
of which the best five were inspected. The best model obtained from 
CYANA for the finally accepted NOE input was subjected to structural 
verification by molecular dynamics (MD) simulations performed with the 
AMBERll software package. ^^^^ The DNA molecule was parameterised 
with AMBER-fflO, the Ru" ion coordination sphere was obtained with 
an earlier described procedure,^^^'^^^ and the AMBER GAFF force 
field^^^^ was used for the parameterisation of the remaining atoms of A,A- 

B. By using standard protocols, the complex of DNA with A,A-B was sol- 
vated by explicit solvent water (TIP3P)^^^^ and neutralised by sodium 
ions, energy minimised, heated and equihbrated during 2 ns with a series 
of pull-and-relax cycles, where NOE restraints were switched on and off, 
by means of steered MD.^^^ "^^^ Both this equilibration stage and the fol- 
lowing productive runs were performed by using a replica-exchange MD 
approach (REMD)^^^^ with a temperature step of 1° covering a tempera- 
ture span of 300-309 K, providing an ensemble of ten paraUel MD runs. 
The productive MD simulations used for the analysis consisted of ten tra- 
jectories of 5 ns each. These were started from the equilibrated struc- 
tures, and simulated without any restraints from NOEs or for hydrogen 
bonds, by using constant pressure and temperature. More details on the 
MD simulations are provided in the Supporting Information. 

Data analysis : Distances corresponding to observed NOEs were averaged 
over MD runs as given in Equation (1): 



d.o.={^j:df^ (1) 

where i is a snapshot of the MD trajectory, A'^ is the number of snapshots 
in this trajectory and di is the distance in this snapshot between two pro- 
tons of the NOE. The interaction energy between A,A-B and the DNA 
was approximated by Coulomb potentials, by using a relative permittivity 
of 20^^^^ between the two Ru ions and the phosphates, Lennard- Jones po- 
tentials between aU intermolecular pairs of atoms and hydrogen bonds 
between donor hydrogen atoms and acceptor atoms closer than 3 A. For 
the former, point charges of +2 centred on the Ru ions and of -1 cen- 
tred on the phosphorus atoms were assumed. The Lennard- Jones param- 
eters were taken from AutoDock 3.0.5,^"^^^ by using the same parameters 
for the phosphorus (missing in the reference) as for sulfur and a cut-off 
of 6 A was used. The large hmit for the hydrogen-bond lengths was 
chosen because the structures investigated represent averaged structures 
over trajectories. Contributions of the intermolecular energy to the sur- 
face created by cutting the DNA (i.e., the extra surface created by ex- 
tracting a 12-mer from a long DNA) consist of interactions from any 
A,A-B atom to the base atoms of the top base pair. For every such inter- 
action, the angle of the interaction direction and the plane of the base 
were calculated: if this angle is larger than 45°, then the interaction is 
considered to involve the surface created by cutting the DNA (i.e. , an in- 
teraction not possible in the presence of a long DNA). 
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